New harmonic number identities with apphcations 



Roberto Tauraso 
Dipartimento di Matematica 
Universita di Roma "Tor Vergata" , Italy 
Q ■ taurasoSmat . uniroma2 . it 

o ■ 

(N 



http : //www.mat .uniroma2 . it/~tauraso 



O 
D 

Q I Abstract 

CO ' We determine the explicit formulas for the sum of products of homogeneous multiple 

harmonic sums J2k=i TVj=i -^fc({l}'^^) when J2]=i — 5- We apply these identities to 
the study of two congruences modulo a power of a prime. 



1 Introduction 

Let s = {si, S2, ■ ■ ■ , Sd) he a vector whose entries are positive integers then we define the 
multiple harmonic sum (MHS for short) for n > as 

■ l<fci<fc2< -<fcd<n ^ ^ 

t:J- ■ We call d and |s| = J2i=i its depth and its weight respectively. This kind of sums has a 

I long history, see for example [1], [12] and their references. 

CN ' In this note we present an algorithmic procedure to determine a closed formula for 

fc=ii=i 

^ ■ involving products of MHS evaluated at n and of total weight less or equal to Yll=i 

H ; 

(i) The first step is to expand recursively any product Hn{s) ■ Hn{t) as a non-negative 
integral combination of MHS of weight |s| + |t|: 

Hn{s)-Hn{t) = 

r6s L±J t 

E Hn{si,r)+ Hnivi,r)+ ^ F„(si + ti,r) 

res' L±J t res L±J t' res' L±J t' 

where s l±j t is the so-called stufRe product (or quasi-shufRe product) of s = (si,s') 
and t = (ii, t') (note that if t is empty then s l±j t = s and Xlres l±j t simply 
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(ii) The second step is to sum over k any MHS by using recursively the following rule: 
fc=i fc=ij=i •' j=i k=j j=i 



(n + l)F(s) 



n 

[, k=l 



where s = {s',Sd)- It is easy to see that the MHS involved in the final formula have 
weight less or equal to |s|. 

(ill) The last step is to simplify the formula by collecting terms by using the stuffle product 
relations. This can be done in various ways and the choices depend on the application 
of the identity. 

This procedure will be employed in the next section and it will generate a large bunch of 
identities. We give two applications of these identities (we hope that many more will come 
out in the future). The first is about the sum 



E(-ir"; 



k=0 

which seems to have a closed form only for a very few values of the parameter a € Z, namely 
for —1 < a < 3. In [2], Cai and Granville studied a related congruence by showing that for 
any prime p > 5 



fe=0 

Here we present the following extension (see also and [3] for similar results): 
Theorem 1.1. Let p > 5 be a prime then for any a ^"L 
p-i 

k=o ^ ~ ^ ap-1 \ 6 

where Xp = — o-n-d Bn denotes the n-th Bernoulli number. 

Moreover, by using the previous theorem for a = — 2 we get 
Corollary 1.2. Let p > 5 be a prime then 

p-i 



2Zl[u)=-^^''^ (modp). 

k=l ^ ' 



which improves the (mod p^) result contained in [9]. 
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2 MHS: identities 

By following the procedure introduced in the previous section, we found an explicit formula 
for any sum of products of homogeneous MHS like Hk{{lY) up to order 5 ({1}"' means that 
the number 1 is repeated d times): 



J2Hk{l) = {n+l)Hn{l)-n, 

k=l 
n 

k=l 

n 

Hl{l) = (n + l)Hl{l) - (2n + 1)F„(1) + 2n , 

k=l 

Y Hk{{lf) = {n + l)H^{{lf) + n {h^{1) - \hI{1)^ + ^i/„(2) - n , 

Y Hkmk{{lf) = {n + + (3n + 1) (k^{1) - \hI{1)) + '^H^{2) - 3n . 
k=l ^ ^ 

J2 Hl{l) = (n + l)Hl{l) + (6n + 3) - \hI{1)\ + \h.^{2) - 6n. 

A;=l ^ ^ 

It's interesting to note that the formulas for X]fe=i ^ki^) when r = 1, 2, 3 appear as Entry 8 
at page 94 in [T]. To illustrate the procedure we show how to obtain Yyk=i -f^fc(l)-f^fc({l}^)- 
By (i) 

= 3i/fc({l}3) + Hki2, 1) + Hk{l, 2). 

By (ii) 



k=l 

n 

Y Hk{2, 1) = (n + 1)F„(2, 1) - nHn{2) + Hn{l) , 

k=l 
n 

Y Hk{l, 2) = (n + l)Hn{l, 2) - Hn{{lf). 

k=l 

Finally by (iii), since by (i) 

Hn{l)Hn{{lV) = 3i/„({l}3) + Hn{2, 1) + 2) and 2F„({1}2) = hI{1) - F„(2), 

we get the formula given above. 

The formulas when the total weight is 4 are contained in the next table: the sum 



Yfk-{n+ l)fn, 



k=l 



where /„ is an entry in the first row, is equal to the sum of the entries of the first column 
each multiplied by the linear polynomial an + b contained in the intersection of the chosen 
row and column. 
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l^k=l Jk- [n+ L)Jn 


IT (f-i \4\ 












— Tl 




/In 1 








—n 


-2n - 2 


-2n-l 


-2n - 3 


-4 




—n 


1 


—n — 1 


1 


3 


\Hn{l)Hn{2) 


n 


2n 


2n + 1 


2n + l 










1 





1 


2 


n 


1 


6 


4 


12 


24 



The large table which gives the formulas when the total weight is 5 is in the Appendix. 



3 MHS: congruences 

Among the various known results about MHS modulo power of a prime, the following ones 
will be crucial for us: for any prime p > 5 

Hp^iil) = 2p^Xp (mod /) , 

Hp^i{2) = -4pXp (mod p3) , 

Hp^i{3) = (mod p^) , 

Hp.i{l,2) = -6Xp (modp2), 

i7p„i(4) = i7p_i({l}2, 2) = Hp^iil, 3) = (mod p). 

where Xp = — ™d i?„ denotes the n-th Bernoulli number (see [7] for the MHS 

of depth 1, see [12] and [3] for all the MHS of depth > 1 with the exception of Hp^i{l,2) 
modulo p'^ which has been established in |10j). 

Note that every homogeneous MHS can be expressed in terms of MHS of depth 1. More 
precisely (see Theorem 2.3 in [3]): given a positive integer d then for any unordered partition 
A = {Ai, A2, . . . , A,.} of d there is some integer cx such that 

r 

d\Hn{{l}'')= Cxl[Hn{Xi). 

For example: 

2Hn{{lV) = Hlil) - H^{2) , 

6H^{{1}^) = Hl{l) - 3Hn{l)Hn{2) + 2Hn{3) , 

24//„({l}4) = H^{1) - ml{l)Hn{2) + 8Hn{l)Hni3) + 3,Hl{2) - 6F„(4). 
Hence, for any prime p > 5 

Hp_i{{lY) = 2pXp (mod p^) , 
Hp^i{{lf) = (mod p2) , 

Hp^i{{l}^) = (mod p). 

Therefore, by the previous identities and congruences, it follows that for any prime p > 5: 
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Efc=igfc(l) = -(p-l)+ V^P (mod p5 

Elll Hkiim ^ (P - 1) + (4 - 2p)p^Xp (mod p^ 

Hl{l) = 2{p - 1) + (2 - 4pyXp (mod p^ 

Ed Hk{{l?) = -(p - 1) + (2 - Ap)pXp (mod p3 

Yl=\ Hk{l)Hk{{lY) = -3(p - 1) + {-Qp)pXp (mod 

Efc^l ^l(l) = -6(P - 1) + (-2 - Qp)pXp (mod p3 

E^;l Hum') = (P - 1) + (-2p)Xp (mod p2 

EU ^ 6(p - 1) + (-6)X, (mod p2 

E^;} iffc(l)i?fe({l}3) ^ 4(p - 1) + (-2p)Xp (mod p2 

//|(l)//fc({l}2) ^ 12(p - 1) + (-6 + 2p)X, (mod p2 
Efc=l gfc(l) = 24(p - 1) + (-12 + 8p)Xp (mod p^ 

EU^fc({l}')^l (modp) 

EU^fc(l)^fc({l}')=5 (modp) 

EU = 10 + 6Xp (mod p) 

E^l} if|(l)/7fc({l}3) ^ 20 + QXp (mod p) 

Ed i^fc(l)i/|({l}2) = 30 + 18Xp (mod p) 

Hl{l)Hk{{lY) = m + 30Xp (mod p) 

Yl~=i Hl{l) = 120 + 60Xp (mod p) 



Note that Efc=i -^fc(l) (mod for r = 1, 2, 3 have been estabhshed by Z. W. Sun m [8]. 

4 Proof of Theorem 1.1 and Corollary 1.2 

Proof of Theorem 1.1. Assume that p > 5 is a prime, then for A; = l,...,p — Iwe have that 

h 5 

{-^f{^~^) = t{U--)^i + iZ^-pym{m (mod/). 

Hence 

T-i+E(-p)^e(:) e (a, ' ,)ri^^(m^o (mod/) 

where P(j, r) is the set of the integer partitions A of j into r parts. 
By summing over k we find 

fe=0 ^ ^ i=l r-=l ^ \evu,r) ^ k=li=l 

Finally, by the congruences estabhshed in the previous section we can compute 

p— 1 r 

(-PrEn^'^({l}'') (mod/) 
k=l 1=1 
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for any partition A G 'P{j, r) and we get easily the result. □ 
Proof of Corollary 1.2. In [6] Staver proved that for any integer n>\ 



^ 1 f2k\ _ f2n\ 2n + 1 ^ /n - 1^ 
fc=l ^ ^ ^ ^ fc=0 ^ 



Letting n = p, by Theorem 1.1 for a = — 2 we have 

^"^ 1 (2k\ 1 /2p\ / 2p + 1 ^ /p - 1^ 



2^2p-lV^^ 8,^\ \_ 16,^ ._^.4^ 



where in the last step we used the fact that i^^Zi) = 1 (mod p^) by Wolstenholme theorem. 

□ 
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Appendix 



ELlA-(" + l)/n 


Hum'') 




i^„(l,l)^n({l}') 


Hlil)Hni{lf) 




n 


5n + l 


lOn + 3 


20n + 7 


\Hn{2) 


n 


3n+ 1 


4n + 3 


6n + 5 




n 


2n+ 1 


n 


2n + 1 


\Hn{l)Hn{2) 


—n 


-3n - 1 


-An - 1 


— 6n — 3 


Hn{l,2) 








-1 


-1 




n 


n+1 


-1 


-1 




—n 


-{n+l) 


-(2n- 1) 


1 




n 


3n+l 


4n + 1 


6n + 3 


lHn{l)Hn{3) 


—n 


-2n- 1 


— n 


-2n - 1 


^n(l,3) 






















1 


1 


n 


-1 


-5 


-10 


-20 



ELl /fc - + l)/n 




Himni{i}') 






30n + 12 


60n + 27 


120n + 60 


\Hn{2) 


6n + 8 


6n + 13 


20 




-3 


-6 


-15 


\Hn{l)Hn{2) 


-6n - 2 


— 6n — 3 





Hn{l,2) 


-3 


-5 


-10 


i^n(4) 


-2 


-3 


-4 


lHU2) 


-2n + 4 


9 


20 


\Hl{l)Hn{2) 


6n + 2 


6n + 3 


















1 


2 


5 


Hnm\2) 


3 


5 


10 


n 


-30 


-60 


-120 
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